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ABSTRACT 


Two triglycine sulphate crystals were grown from an aqueous solution in 

Spacelab 3 aboard a Space Shuttle. Using a diffusion coefficient of 2 X 10 ^ 

2 

cm /s, a computer simulation gave reasonable agreement between experimental 

and theoretical crystal sizes and interferometric lines in the solution near 

the growing crystal. This diffusion coefficient is larger than most measured 

-3 

values, possibly due to fluctuating accelerations on the order of 10 g 
(Earth’s gravity). The average acceleration was estimated to be less than 
10 ^ g. At this level buoyancy-driven convection is predicted to add 
approximately 20% to the steady-state growth rate. Only very slight 
distortion of the interferometric lines was observed at the end of a 33 hour 
run. It is suggested that the time to reach steady state convective transport 
may be inversely proportional to g at low g, so that the full effect of 
convection was not realized in these experiments. 

* Present address: Samsung Advanced Institute of Technology, Kihung-Eup, 

Kyungki-Do, South Korea. 
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1. Introduction 


In May of 1985 a layer of triglycine sulphate (TGS) was grown from 
aqueous solutions onto two seed crystals in Spacelab 3 aboard a Space Shuttle. 

The objectives of the experiment were to determine how crystal growth from 
solution progresses in a nearly convection-free environment, the influence of 
low g growth on crystal properties, and the extent of buoyancy-driven 
convection in the small fluctuating accelerations that are present in the 
Shuttle . 

In the absence of convection, crystal growth from solutions will rapidly 

become slower as the adjacent solution is depleted of solute, unless the 

growth temperature is lowered to compensate. However if the temperature is 

lowered too fast the growth rate will reach the level at which solvent 

inclusions are formed (1). In order to estimate the optimal temperature 

program prior to the Spacelab 3 mission, a one dimensional spherical computer 

-5 2 

model was developed and used with a diffusion coefficient D of 5 X 10 cm / s 

_ i 

in the solution (2,3). Later the computations were repeated with D * 1 X 10 
2 

cm /s using both the spherical model and a newly developed two-dimensional 
cylindrical model (4). It was estimated that only about 10% of the 
supersaturation of the bulk solution is required to drive the interface 
kinetics; i.e. the growth is 90% diffusion-controlled. 

Here we report on a comparison of the computations using the cylindrical 
model and the experiments. Estimates are made of the buoyancy-driven 
convection that occurred and its development with time. 

PRECEDING PAGE BLANK NOT FILMED 


291 


2 . Spacelab 3 experiments 


The growth experiments were performed in the Fluid Experiment System 
(FES) designed and constructed by TRW Systems in Redondo Beach, California, 
with final debugging by NASA's Marshall Space Flight Center (MSFC) in 
Huntsville, Alabama. Based on the original concept of Lai, the growth cell in 
the FES was designed so that a seed crystal is cemented onto a controlled 
temperature rod (sting) (5,6). An optical system provided real time schlieren 
video, transmission holograms and reflection holograms of the entire interior 
of the cell (7,8). 

Quantum Technology of Sanford, Florida provided 3.4 mm thick (001) 
cylindrical seeds for the two growth runs. For Cell 2 the seed was 15 mm in 
diameter and for Cell 3 10 mm diameter. In each cell the seed was covered by 
a polymer diaphragm cap to prevent the solution from contacting the seed 
prematurely. The cell was filled with a TGS solution of concentration 
equivalent to saturation at 45 C (9). During storage of the cell at room 
temperature TGS crystallized throughout the cell. In Spacelab 3 the cell was 
slowly heated to 70 C and held at that temperature for several hours to 
redissolve this material. A pump provided circulation. The temperature was 
reduced to 46.5 C for Cell 2 and 46 C for Cell 3, the circulating pump turned 
off, and the cap withdrawn from the seed crystal. The cap in Cell 2 had 
leaked and spurious crystals were seen upon cap removal. Several of these 
later grew around the periphery of the crystal, while one floated in the 
solution. 

At cap withdrawal the solution temperature was above its saturation 
temperature so that dissolution was expected. The schlieren video confirmed 
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that dissolution was taking place as planned. After 30 minutes for Cell 2 and 
40.8 minutes for Cell 3, the sting temperature and the cell wall temperature 
were slowly reduced. Dissolution continued for some time before the schlieren 
showed the beginning of growth. Figures 1 to 3 show the temperature vs. time 
data. Note that during the Cell 2 experiment the ground crew called for a 
temperature increase because they feared the growth rate was becoming too 
large . 

During growth, accelerometers in FES measured accelerations as high as 

-3 

0.1 of Earth* s gravity (g), with average fluctuations of 10 g and a maximum 

in the power spectrum at 17 Hz. The average DC acceleration level was 
-6 

approximately 10 g. 

After growth the cap was replaced over the crystal. The sting was 
withdrawn from the cell and excess solution blotted off with soft tissue. 

Later the crystal was pried from the sting, measured and photographed (Figures 
4 and 5). While crystal properties will be reported in detail elsewhere (27), 
it is noteworthy that the boundaries between the seed crystals and the 
space-grown layers could not be seen using optical microscopy. (Earth-grown 
crystals typically form inclusions at this boundary, with dislocations 
emerging from the inclusions into material grown on top of them. (e.g. Ref. 
10)) Presumably this is a reflection of the very slow transition from 
dissolution to growth that occurs under diffusion-controlled conditions. The 
top surface of both crystals was concave, indicating more rapid vertical 
growth around the periphery. Around the edge (001) and (010) facets formed. 
From one edge of the Cell 2 crystal a thin transparent sheet had grown over 
almost the entire surface. Such a feature has never been observed on Earth, 
possibly because such a thin sheet would be bent or broken due to gravity. 
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Holograms were developed and converted into views of the growth cell and 
interf erograms of the solution, both parallel and perpendicular to the sting. 
Typical interf erograms of Cell 3 are given in Figures 6 to 8 . The 
interf erograms from Cell 2 are badly distorted because of the perturbation of 
the concentration field by growth of the spurious crystals. If cylindrical 
symmetry is assumed the interf erograms can be converted to refractive index 
maps (and vice versa), permitting the comparison of theory and experiment 
described later. 

The side views of Cell 2 permitted R. Naumann at MSFC to measure the size 
and movement rate of the floating spurious crystal. Use of Stokes law with 
the viscosity of the solution and the densities of solution and crystal (9) 
yielded an acceleration level of 4 X 10 ^ g (8). This must be regarded as a 
lower limit because movement of the crystal parallel to the line of view was 
not determined. Furthermore the solution was considered to be stationary, 
whereas some buoyancy-driven convection may have been present. Stokes law 
assumes that the solid object is spherical and that the solution is of uniform 
density, while the actual crystal was not spherical and the solution near the 
crystal had a lower density because it was depleted in solute due to growth. 


3. Mathematical model 


Details of the mathematical model, including all of the FORTRAN computer 
programs, are given in Reference (4). Figure 9 shows the geometry used for 
the numerical model of diffusion and heat conduction during growth. Heat and 
mass transfer due to buoyancy-driven convection were ignored. (The 
consequences of this are discussed later.) 
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Cylindrical symmetry was assumed. While the sting and the seed are 
cylindrical, the cell walls constitute a 10 X 10 X 10 cm cube. In addition a 
temperature sensor runs parallel to the sting and a vacuum line attaches to 
the sting, as can be seen in Figure 8. Nevertheless a cylindrical model is a 
reasonable approximation because the refractive index field in the solution is 
determined almost entirely by the concentration field, and this varies 
markedly only near the crystal. 


Because the partial molar volume of TGS is different in the crystal and 
in the solution, the growth generates a gentle flow of solution toward the 
crystal. This flow carries both heat and mass. However examination of the 
terms in the transport differential equations revealed that this flow 
contributes only about 0.2% to the mass transfer rate and 0.001% to the heat 
transfer. Consequently these contributions were ignored. If we assume 
constant properties in the solution this leads to the differential equation 
for mass transfer: 



8 , 8Cv 

8? <r 9?' 



8C 

8t 


and for heat transfer: 


e 


a , aiv 
37 (r 97 ) 
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8T 
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( 2 ) 


where the symbols are defined in the Table of Nomenclature. The temperatures 
at the wall and sting tip were taken to be prescribed functions of time 
(Figures i to 3). Heat flux balances and continuity of temperature were used 
at each phase boundary. The normal concentration gradient at sting, cap and 
wall surfaces was set equal to zero because of the absence of diffusion 
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through these boundaries. Along the cylindrical axis (r - 0) the radial 
concentration and temperature gradients were set to zero because of symmetry. 
L'Hopital's theorem was applied to eliminate the singularity in the convective 
terms at r ■ 0. 

The heat transfer and mass transfer are independent of one another except 
at the surface of the crystal, where they are coupled by the relations for 
interface kinetics : 

v c = K(T e - Ti) (3) 

and solubility: 

« 0.237 exp( - 1679/T e )mol/cm^ (4) 


The interface kinetics constant K was derived from the data of References 

(10,11), i.e. for a (001) face a growth rate of 1 mm/day for an interfacial 

undercooling of 1 K and 2.4 mm/day for a (010) face. Since the top face in 

the flight experiments was (001), K = lmm/day.K was used for that face. For 

the side surface and for dissolution on the top surface equilibrium was 

assumed, i.e. infinite K so that T. « T . 

l e 


The relationship between the concentration gradient at the crystal 
surface and the growth rate were found from a material balance taking into 
account the Stefan or crystallization flow (12) and the absence of solvent 
(water) in the crystal: 


D(VC)i 

v c 55 7 =7 

c CcU-CjV) 


(5) 
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The crystallization flow gives rise to the term C^V, which equals the 
volume fraction of solute at the crystal surface (approximately 0.2 in these 
experiments). It is negligible only if the solubility of a salt is 
negligible, which is most definitely not true for TGS. 

The foregoing set of equations were solved using the Multiple Alternating 
Direction Implicit (MADI) finite difference method (13). It is a modification 
of the Alternating Direction Implicit (ADI) method in which new values of 
concentration (or temperature) are calculated for each row by implicit methods 
using values in adjacent rows from the preceding time step (14). In the next 
time step of an ADI computation, new values in each column are calculated 
using values in adjacent columns from the preceding time step. In MADI the 
new value at a point is taken as the average of values calculated from 
adjacent rows and from adjacent columns. MADI allows larger time increments 
to be used than ADI with the same accuracy, thereby reducing the computation 
time required. 

A time increment of 10 s was used for the temperature calculations and 20 
s for the concentration calculations. For isothermal dissolution a negligible 
difference was found for a time step ranging from 0.5 s to 30 s for 
concentration computations. 

Because the concentration gradient was expected to be largest near the 
crystal, a finer mesh was used there as shown in Figure 10. Because of the 
assumed axisyrametry, only one-half of the cell needed to be solved. 

During growth the crystal surface moves relative to the cell, so that it 
is generally between mesh points. In order to track the top crystal surface 
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and the conditions there, the following operations were performed at each time 

step. First the interfacial temperature was found by interpolating between 

temperatures at adjacent mesh points in the crystal and the solution. This 

left three unknowns in Equations 3, 4 and 5; C., V and T (assuming that the 

i c e 

interfacial concentration gradient can readily be calculated once is 
known). After eliminating by setting Equations 3 and 5 equal to one 
another and substituting for from Equation 4, the combination was solved 
for by successive approximations. Knowing was calculated from 

Equation 4 and from Equation 3. 

For the top surface during dissolution and for the side crystal surface 

during both dissolution and growth K was taken as infinite, and so T - T. C, 

e 1,1 

was calculated from Equation 4, and from Equation 5. 

In the computation of concentration and temperature fields the crystal 
was maintained in a cylindrical shape by determining new interfacial positions 
for the entire top surface using the growth rate computed for the midpoint 
between the center and the periphery. For the side surface the growth rate 
halfway down was used. This procedure amounts to ignoring the effect of the 
deviation from cylindrical shape on the heat and mass transfer, resulting in a 
slight underestimation of the growth rate at the periphery (top corner). In 
comparing theory and experiment the crystal shape vs. time was computed from 
the calculated growth rate at each mesh point. 

4 . Comparison of experiment and theory 

Initial comparison of theory and experiment were made for Cell 3 because 
of the problems with the Cell 2 experiment (cap leakage during storage, 
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spurious crystals, and crystal sheet grown above top surface of seed). Figure 

11 shows a comparison of the size of the actual crystal received on Earth with 

“5 2 

computed sizes for D * 1, 2 and 3 X 10 era /s. The best fit with crystal 

-5 2 

thickness (height) was for D * 2 X 10 cm /s, and so this value was used for 
subsequent computations. Note that theory and experiment both produced a 
concave top surface, i.e. the growth rate was higher at the periphery. 

Likewise the growth on the side of the crystal was several times larger than 
on the center of the top face. 

Figure 12 shows isoconcentration lines for Cell 3 at the end of the 40.83 
minute isothermal dissolution period. The solute concentration is highest at 
the crystal surface. As the temperatures of the wall and the sting are slowly 
decreased, dissolution continues for a considerable time, albeit at a 
continually declining rate. Eventually growth begins and the solute 
concentration at the crystal surface declines. This leads to a maximum in 
concentration some distance away from the crystal in the solution, as shown in 
Figure 13. Figure 14 shows the iso-concentration lines at the end of the 
run. Note that there is no significant change in concentration at the cap, 
the cell walls or even the side of the sting. Thus the original assumption of 
cylindrical symmetry is reasonable. 

As the temperature of the wall and the sting tip are reduced after the 
initial isothermal dissolution period, heat flows from the solution to these 
surfaces, resulting in a maximum of temperature in the center of the cell as 
shown in Figure 15. This maximum persists for some time after the end of 
ramping down the temperature, until eventually the solution is again 
isothermal . 
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Using the data of Reference (9) we computed refractive index maps from 
the temperature and concentration data. The variation in concentration has a 
much larger effect on refractive index than does the variation in temperature. 

Thus the refractive index lines look much like the isoconcentration lines. 
Figure 16 shows a refractive index map at the end of the isothermal 
dissolution, where only concentration varies in the cell. Figure 17 shows 
refractive index during the programmed cooldown period, at a time showing a 
maximum in concentration and temperature in the middle of the cell. For long 
times where the solution is again isothermal, the refractive index maps look 
exactly like the concentration maps. 


Interf erograms wore calculated from the refractive index field in order 
to compare with those obtained from the Spacelab 3 holograms (see Figures 6 to 
8). Since the laser beam came in the side of the cell parallel to the top 


ft 


crystal surface, 
direction (y) at 


it was necessary to compute the optical pathlength in that 

A 


every position using: 


l = /n dy (6) 

Bright interference lines pass through points at which the optical pathlength 
has changed by an integral number of wavelengths of the helium-neon laser 
(0.633 microns). Dark lines occur between the bright lines at half- integral 
wavelength changes in optical pathlength. Results are given in Figures 18 and 
19 (for clarity, every other fringe is deleted in Figure 19). Comparison of 
the experimental results in Figure 6 with the calculated results in Figure 18 
for the same time give perfect agreement in the number of fringes above the 
top crystal surface and the number intersecting the top crystal surface. At 
the time represented in Figures 7 and 19, the number of fringes along the top 
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surface of the crystal is the same from experiment and computation. Moving 
straight sideways from the tcp surface the number of fringes from the edge of 
the crystal to a point directly above the edge of the sting insulator is 48 
from the computations and 46 from experiment. 

Finally Figure 20 shows computed and experimental crystal shapes 
resulting from the Cell 2 experiment. In this case the amount of growth on 
the side of the crystal is about the same from experiment and computation, 
while the computed growth on the top face is over twice as much. In view of 
the crystal sheet that grew over the top surface, less growth is not 
surprising because the sheet would have prevented diffusion from the bulk of 
the solution. In addition there is considerable uncertainty in the initial 
position of the crystal because of the cap leakage. Because of the presence 
of spurious crystals during cap retraction it is probable that solution had 
leaked in about the time the cell was filled, causing crystals to nucleate and 
grow during the long storage at room temperature. The seed crystal must have 
grown a little thicker at the same time. During the 70 C dissolution period 
in Spacelab 3 both the spurious crystals and the seed crystal would have 
dissolved, with the total dissolution exceeding the amount that had 
crystallized during storage. Upon cooling to 46.5 C regrowth and possibly 
nucleation would have occurred. 


5. Convection 


While the accelerations in orbiting spacecraft are small, they are not 
zero, as indicated in Section 2 above. Because the density of the solution 
was not constant during the Spacelab 3 FES experiments, some buoyancy-driven 
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convection was to be expected. We need to separate the effects of 
acceleration into those due to fluctuations and those arising from the much 
smaller steady background. We first consider fluctuating effects and compare 
expectations with experimental results. 

As we show shortly, the time constant for fluid mechanics is on the order 
of 1 to 100 minutes in an FES cell. Since the fluctuating components of 
acceleration had much shorter periods than this, one would not expect any 
direct evidence of convection such as rapid shifts in the schlieren pattern. 
Indeed rapid shifts in schlieren were not observed and observable shifts in 
interf erograms were not detected even after 1 day of growth. (Double exposure 
holograms did sometimes reveal local changes in a few seconds.) It is 
conceivable that a fluctuating acceleration could cause minute movements of 
solution sufficient to influence mass transfer, which is very sensitive to 
convection for liquids. Evidence for this might be found in the diffusion 
coefficient required to bring experiment and convection-free theory into 
agreement. If such an effective diffusion coefficient exceeds independently 
measured values, we might conclude that convection affected the mass 
transfer. 


As noted earlier, the effective diffusion coefficient found for these 
-5 2 

experiments is 2 X 10 cm /s. This value equals the average found by Reddy 

using a porous diaphragm technique for temperatures ranging from 25 to 45 C 

and concentrations up to 1.2 molar (15). However the scatter in the data was 

on the order of + /- 100%. Using a similar technique Kroes and Reiss obtained 

-5 2 

for very dilute solutions a value of 1.3 X 10 cm /s at 45 C (9). Using a 
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ref ractoraetric technique Novotny, Moravec and Sole obtained 0.7 X 10 cm /s 
for a saturated TGS solution at 45 C. While the foregoing indicates that 
fluctuating acceleration may increase mass transfer rates, the evidence is 
weak. A computer simulation would be most helpful, but a supercomputer is 
needed for a full three-dimensional time-dependent coupled fluid mechanics / 
mass transfer solution. 

-6 

The steady background acceleration, on the order of 10 g, is expected 
to generate convection after some period of time (as discussed later). We can 
use existing correlations (16-21) to estimate the steady state contribution of 
this convection to the growth rate. (These correlations are for 
buoyancy-driven natural convection about spheres (16), cylinders (17, 21) and 
vertical flat plates f 18-20) with low Grashof number.) The ratio of the 
contribution of convection to the growth rate to that for diffusion is given 
by : 

v conv Sh CO nv AC ( 7 ) 

v dif f = d(VC)i 

where Sh is the contribution of convection to the Sherwood number from a 

conv 

correlation. The concentration difference &C between solution at the crystal 

3 

surface and far from it is taken to be 0.057 mmol/cm , the crystal diameter 

4 

1.5 cm, and the interfacial concentration gradient 0.44 mmol/cm . Using the 
properties given in the Table of Nomenclature for TGS solutions we calculate a 
value of 0.14 for the Grashof number Gr at 10 ^ g and 650 for the Schmidt 
number Sc. From Equation (7) and the correlations we estimate a convective 
contribution to the growth rate ranging from 13 to 22% of the diffusive 
contribution. 


1 
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The steady state fluid velocity due to convection is on the order of 


( 22 ): 


V = v Gr/d 


for a Grashof number much less than unity, and 
V = vGr^/d 


( 8 ) 


for a Grashof number much greater than unity. With the values given above (Gr 
* 0.14) we use Equation (8) to estimate V * 4.5 cm/hr, which is sufficiently 
large that it should have been obvious in both the schlieren pattern and the 
intef erograms for Cell 3. The fact that distortion was not observed until 
very late in this run indicates that steady state convection was not obtained. 

We must ask ourselves how long it takes to reach steady state natural 
convection at low acceleration levels. 


Three-dimensional time-dependent convection computations were carried out 

using the discrete element method (23). Unfortunately computer requirements 

were extremely heavy so only short periods of real time could be simulated. 

In order to isolate initial transient diffusion effects from transient 

convection effects, diffusion was allowed to proceed for 46 hours in the 

absence of convection for a concentration difference of 0.057 mmol/cc between 

the top crystal surface and the bulk solution. Figure 21 shows the maximum 

-4 

solution velocity in Cell 2 after turning on an acceleration of 10 g (Gr = 
14). If the velocity continued to increase at the same rate, it would take 
0.36 hours to reach the steady state velocity of 0.032 cm/s estimated from 
Equation (9). However in practice the rate of velocity increase with time 
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would diminish, so that more than 0.36 hours would be required to reach steady 
state. The question is, how much time would be required at 10 ^ g? If the 
concentration field is not influenced by convection, which is not likely to be 
true for large Sc, the hydrodynamic transient time is independent of g level. 
However for the coupled mass transfer and fluid mechanics that occurs in 
natural convection, this may not be true. 

In the literature we located three boundary layer theoretical treatments 
(good for large Gr) of transient natural convection from a vertical flat plate 
(24-26). All of these treatments agree on the Gr dependence of the transient 
time, but only one (24) examined the Sc dependence. The time required to 
reach steady state after turning on a temperature or concentration difference 
is approximately (24): 

t = 5.24(0.95 + Sc ) 1 / 2 Gr 1 / 2 x 2 /v 


from which we calculate 1.69 hours for the Cell 2 conditions examined above 

-4 -6 

with an acceleration of 10 g and 16.9 hours at 10 g. However for low Gr, 
as we have for 10 ® g, it is likely that the 1/g dependence of transient time 
is greater than the one-half power. Since the g dependence of velocity goes 
from 1/2 power to 1 as Gr becomes small (Equations 8 and 9), it is likely that 
the same thing happens for transient time. If we assume a time dependence of 
1/g, then we estimate a transient time of 169 hours to achieve steady state 
conditions. This is far longer than the length of the Spacelab 3 experiment! 
But is this reasonable? 
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6. Transients 


In order to gain some theoretical guidance on the various transient 
times, we turn to the governing partial differential equations. The 
Navier-Stokes equation for hydrodynamics with the usual Boussinesq 
approximation is: 


+ V • VV = v V2y - B gAC 
d t 


(id 


The mass transfer equation with convection is: 
12 = DV 2 c - V • VC 


( 12 ) 


The convective heat transfer equation for low velocity flows is: 


pc 


9T 

p at 


kV 2 T - pCpV • VT 


( 13 ) 


We non-dimensionalize these PDE by dividing all lengths by a 
characteristic length L, all concentrations by a characteristic concentration 
difference C, all temperatures by a characteristic temperature difference 
^T, all velocities by a reference velocity V , and all times by a reference 
time t^. Thus the above equations become: 


L 3 V r av L 2v r 2 

+ — V • vv 


LV 


v‘t. 


ax 


- V2 v - Gr 


( 14 ) 


where v - V/V^ is the dimensionless velocity. The first term in this equation 
is for acceleration, the second for inertia, the third for viscous effects, 
and the third comes from buoyancy and drives the convection. Similarly the 


306 



dimensionless mass transfer equation is: 


L -— M = v2<(, 

t r D 


V r L 


v<{> 


( 15 ) 


The first term is the transient terra, the second the diffusive transport term 
and the third the convective transport term. The dimensionless heat transfer 
equation is. 


l 2 P c P 94) 
t r k 9 t 


V 2 4) - 


V r Lpc p 

k 


v • Vip 


(16) 


All of the pre-differential expressions and all of the differentials in 
the above equations are dimensionless. The non-dimens ional ising parameters 
must be chosen such that every differential in the PDEs is of order one. Thus 
a logical choice for L is the crystal diameter, for & C the difference in 
solution concentration between the crystal surface and the cell wall, and for 
the difference in temperature between the sting and the cell wall. The 
reference velocity is on the order of the maximum velocity in the system and 
the reference time is on the order of the transient time. The maximum 
velocity and the transient time are what we want to estimate. 

If acceleration is zero, Gr * 0. Noting that viscous decay damps 
convection in the absence of a driving force, we set the acceleration term 
equal to the viscous term to obtain the following estimate for the transient 
time due to fluid mechanics alone: 

= L 2 /v 


(t r)fluid 


(17) 


2 

For a characteristic length of 1 cm and a kinematic viscosity of 0.013 cm /s, 
Equation 17 gives a transient time of 77 seconds. 

In the absence of convection the convective transport terms in the heat 
and mass transfer PDEs are zero and the pre-differential expressions in them 
must be of order one. Thus the transient times for diffusion and conduction 
are : 


<t r > di ff = l2 /° 


(18) 


^r^cond L^pCp/k (19) 

From these we estimate a transient diffusion time of 14 hours and a transient 
conduction time of 19 minutes for our experiments. 

For large Gr, the viscous terra in the momentum Equation 14 is negligible 
compared to the inertial term. Thus the pre-differential expression in the 
intertial term must equal the buoyancy term, from which: 

v r = Gr 1 / 2 v/L (20) 

This is the same as reported in the literature (22). In the mass transfer 
equation 15, convective effects are fully established when the transient term 
is about equal to the convection term. Equating the pre-differentials and 
substituting equation 20 we obtain: 
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( ^comb 


oc 


( 21 ) 



L 2 

v Gr 1 / 2 


1 

1/2 


which is the same functional dependence on g as reported in the literature for 
large Gr (24-26). (Implicit in the foregoing is the assumption that density 


differences arise entirely from concentration differences, which is a good 
approximation in crystal growth from solution.) 


For small Gr, the inertial term in equation 14 is negligible compared to 
the viscous term. Equating the pre-differential expression in the viscous 
term to buoyancy we obtain: 

v r = Gr v/L (22) 


which is the same as given in the literature (22). Substituting this into the 
mass transfer equation and equating transient and convective pre-differentials 
we obtain: 


^ **r ^comb 


_ L2_ 1 

V r vGr “ g 


(23) 


7 . Summary 

Both experiments and theory indicate that the time required to reach 
steady state in buoyancy-driven mass transfer increases as the average 
acceleration level decreases. The same might be true for heat transfer 
controlled free convection, as is common in solidification. Thus some 
experiments performed in space may terminate before the convection reaches 
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steady state levels. Transient computer simulations are needed to confirm 
these predictions . 

The experiments indicate that fluctuating accelerations may significantly 
enhance mass transport rates, even when global convection is not generated. 
Transient three dimensional computer simulations are needed to confirm this 
result . 

The lack of vigorous convection in a space environment makes it easy to 
achieve a gentle transition from dissolution to growth so that solvent 
inclusions are not formed between the seed and the grown layers. Experiments 
on earth have shown that closure of such inclusions tends to result in 
dislocations which propagate through subsequently grown materials and degrade 
properties. It may be that this will permit us to grow crystals from solution 
in space superior to what is possible on earth. More flight experiments are 
needed to adequately test this possibility. 
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Nomenclature 

(Properties used in computations are given in parentheses) 


C 

AC 

C c 

Ci 

c p 


D 

d 

Gr 

g 

&a 

K 

k 


L 

Z 

M 

n 

Sc 

r 

^conv 

T 

T e 


Concentration of TGS in solution (mol/cm 3 ) 

Difference in solution concentration at crystal surface and in 
bulk solution (mol/cm 3 ) 

TGS concentration in crystal (p c /M) (5.23 x 10~ 3 mol/cm 3 ) 

Solution concentration at crystal surface 
Heat capacity - of solution (5 x 10 3 J/kg.K) 

- of crystal (1.5 x 10 3 J/kg.K) 

- of Teflon (1.26 x 10 3 J/kg.K) 

- of Dow Corning 732 RTV adhesive (1.46 x 10 3 J/kg.K) 
TGS diffusion coefficient in solution (2 x 10“ 3 cm 3 /s) 

Diameter of crystal (1 and 1.5 cm) 

Grashof number for mass transfer, g a f3d 3 AC/v 3 
Earth’s gravitational acceleration (9.8 m/s 3 ) 

Acceleration (m/s 3 ) 

Interface kinetics constant (1.16 x 10“® m/s.K for (001) TGS) 
Thermal conductivity - of solution (0.5 J/m.k.s) 

- of crystal (0.68 J/m.k.s) 

- of Teflon (0.41 J/m.k.s) 

- of cement (0.19 J/m.k.s) 

Characteristic length (d or half of cell width) (m) 

Optical path length (m) 

Molecular weight (323 g/mol for TGS) 

Refractive index 

Schmidt number of solution, v/D (650) 

Radial position (m) 

Sherwood number due to convective mass transfer 
Temperature (K) 

Equilibrium temperature corresponding to solubility (K) 
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Temperature at crystal surface (K) 

Temperature difference (K) 

Time since retraction of cap (s) 

Characteristic transient time (s) 

Partial molar volume (181 cm^/mol) 

Velocity (m/s ) 

Crystal growth rate (m/s) 

Reference velocity (m/s) 

V/V r , Dimensionless velocity 
Axial position (m) 

Concentration densif ication coefficient of solution, (-3p/9C)'p/p 
(0.46 cm^/mol) 

Gradient (m“^) 

(C -C^)/AC, Dimensionless concentration 
Kinematic viscosity of solution (0.013 cm 2 /s) 

(T - T^)/AT, Dimensionless temperature 
Density - of solution (1.16 g/ cm^) 

- of crystal (1.69 g/cm^) 

- of Teflon (2.3 g/cm^) 

- of cement (1.03 g/cm^) 
t/t r , Dimensionless time 
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TEMPERA TURE 
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Figure 1. Measured sting temperature vs. time for Cell 3. Wall 
temperature the same within experimental error. 
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Figure 2. Measured sting temperature vs. time for Cell 2. 
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Figure 3. Measured wall temperature vs. time for Cell 2 
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Figure 7. Cell 3 interferogram 16.87 hours after cap withdrawal. 
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Figure 8. Cell 3 interferogram 24.33 hours after cap retraction. Actual 
image of sting and crystal are seen in center. Vertical rod to 
left of sting is thermal sensor. Rod curving down from right 
side of sting is vacuum line. 
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Figure 10. Mesh used for finite difference solution. Numbers 
are distances in cm. 
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Radius (cm) 


Figure 11. Comparison between the final shape of the crystal grown 
in Cell 3: 

(1) Averaged from experimental values in Figure 4. 

-5 2 

(2) Calculated using D - 3 X 10 cm /s. 

-5 2 

(3) Calculated using D * 2 X 10 cm /s. 

-5 2 

(4) Calculated using D * 1 X 10 cm /s. 
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Figure 12. Isoconcentration lines for Cell 3 at end of isothermal 
dissolution period, 40.83 minutes after cap retraction, 

-5 2 

calculated using D » 2 X 10 cm /s. T - 46 C, 

3 

concentration at crystal surface * 1.229 raraol/cra , 

3 

original concentration of solution = 1.209 mraol/cm , 

3 

change in concentration between lines = 0.002 mmol/cra . 
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Figure 13. Computed isoconcentration lines for Cell 3 during linear 
cooldown period, 208.67 minutes after cap retraction. 

3 

Change in concentration between lines = 0.001 ramol/cm . 
Maximum concentration at about Z = 5.4 cm, R ■ 0. 
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Figure 14. Computed isoconcentration lines for Cell 3 at end of run, 
33.97 hours after cap retraction. 

3 

Change in concentration between lines = 0.005 mmol/cm . 
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Figure 15. Isotherms for Cell 3 during linear cooldown period, 330.62 
minutes after cap retraction. Wall temperature and sting 
temperature are 43.9 C. 

Change in temperature between lines * 0.05 C. 






Figure 16. Computed lines of constant refractive index for Cell 3 at 
end of isothermal dissolution, 40.83 minutes after cap 
retraction. Refractive index of solution at top crystal 
surface is 1.38345. 

Change in refractive index between lines = 0.0001. 
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Figure 17. Computed lines of constant refractive index for Cell 3 in 
the linear cooldown period, 208.67 minures after cap 
retraction. Refractive index of solution at top crystal 
surface is 1.38286. 

Change in refractive index between lines « 0.00002. 
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Figure 18. Computed interf erogram for Cell 3 at the end of the 

isothermal dissolution period, 40.83 minutes after cap 
retraction. Compare with experimental interf erogram in 
Figure 6. 
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Figure 19. Computed interf erogram for Cell 3 in the constant 

temperature period, 16.87 hours after cap retraction. 
Every other interference line was deleted for clarity 
Compare with experimental interf erogram in Figure 7. 
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Figure 20. Comparison between experimental and computed crystal 
shape for Cell 2. 

(1) Calculated using D * 2 X 10 5 cm 1 2 3 /s. 

(2) Averaged from experimental measurements. 

(3) Original seed crystal. 
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